Astronomy & Astrophysics manuscript no. 


February 5, 2008 


(DOI: will be inserted by hand later) 





The chemical composition of metal-poor emission-line galaxies 
in the Data Release 3 of the Sloan Digital Sky Survey* 

Y. 1. Izotov\ G. Stasinska^, G. Meynet^, N. G. Guseva\ and T. X. Thuan^ 



Main Astronomical Observatory, Ukrainian National Academy of Sciences, Zabolotnoho 27, Kyiv 03680, Ukraine 
LUTH, Observatoire de Meudon, F-92195 Meudon Cedex, France 
Geneva Observatory, CH-1290 Sauverny, Switzerland 

Astronomy Department, University of Virginia, Charlottesville, VA 22903, USA 



Received ; Accepted 

Abstract. We have re-evaluated empirical expressions for the abundance determination of N, O, Ne, S, CI, Ar and 
Fe taking into account the latest atomic data and constructing an appropriate grid of photoionization models with 
state-of-the art model atmospheres. Using these expressions we have derived heavy element abundances in the ~ 
310 emission-line galaxies from the Data Release 3 of the Sloan Digital Sky Survey (SDSS) with an observed H/3 
flux F(H/3) > 10"^* erg s~^ cm~^ and for which the [O ill] A4363 emission line was detected at least at a 2a 
level, allowing abundance determination by direct methods. The oxygen abundance 12 -I- log 0/H of the SDSS 
galaxies lies in the range from ~ 7.1 {Zq/30) to ~ 8.5 (0.7 Zq). The SDSS sample is merged with a sample 
of 109 blue compact dwarf (BCD) galaxies with high quality spectra, which contains extremely low-metallicity 
objects. We use the merged sample to study the abundance patterns of low-metallicity emission-line galaxies. We 
find that extremely metal-poor galaxies (12 + log 0/H < 7.6, i.e. Z < Zq/12) are rare in the SDSS sample. The 
Q element-to-oxygen abundance ratios do not show any significant trends with oxygen abundance, in agreement 
with previous studies, except for a slight increase of Ne/0 with increasing metallicity, which we interpret as due 
to a moderate depletion of O onto grains in the most metal-rich galaxies. The Fe/0 abundance ratio is smaller 
than the solar value, by up to 1 dex at the high metallicity end. We also find that Fe/0 increases with decreasing 
H/9 equivalent width EW(H/3). We interpret this as a sign of strong depletion onto dust grains, and gradual 
destruction of those grains on a time scale of a few Myr. All the galaxies are found to have log N/O > -1.6, 
implying that they have a different nature than the subsample of high-redshift damped Lya systems with log 
N/O of ~ -2.3 and that their ages are larger than 100 - 300 Myr. We confirm the apparent increase in N/O 
with decreasing EW(H/3), already shown in previous studies, and explain it as the signature of gradual nitrogen 
ejection by massive stars from the most recent starburst. 
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1. Introduction 

Knowledge of the detailed chemical composition of galax- 
ies is fundamental to our understanding of stellar nucle- 
osynthesis and evolution of galaxies. Emission-line galax- 
ies provide an easy way to determine the abundances of 
such elements as He, N, O, Ne, S, CI, Ar and Fe, out to 
redshifts of about 0.4, from an analysis of the radiation 
from their H ii regions in the visible domain. These de- 
terminations are considered more reliable if the electron 
temperatures can be measured directly, using the [O ill] 
A4363/[0 III] A5007 line ratio. Then, the ionic abundances 
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can be derived directly from the strengths of the emission 
lines. The assumptions to derive elemental abundances are 
reduced to a minimum. The weak [O ill] A4363 can only 
be measured in metal-poor emission-line galaxies, where 
the H II regions suffer little cooling and are at high enough 
temperature to produce significant emission of this emis- 
sion line. 

Metal-poor galaxies are the least chemically evolved 
galaxies and thus provide the simplest test beds for theo- 
ries of chemical evolution of galaxies. They also possibly 
constitute a local counterpart to primeval high-redshift 
galaxies. In this spirit, they have be en the subject of 
many studies in the recent past (e.g.. Hzotov fc Thuan L 
l2004bl and references therein). They have also been used 
for the quest of the abun dance of primordial helium (e.g. 
Ilzotov fc Thuan I l2004aL and references therein). Those 
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studies were based on limited samples of very different 
coverages and selection criteria. 

The Sloan Digital Sky Survey (SDSS, lYork et al. L 

1200(1) offers a gigantic data base of galaxies with well- 
defined selection criteria and observed in a homogeneous 
way. In addition, the spectral resolution is much better 
than that of most previous data bases on emission-line 
galaxies. From this data base, it is possible to extract a 
sample of emission-line galaxies with well defined criteria. 

Two papers have already used the SDSS to find ob- 
jects with detected [O lii] A4363. The first one is by 
llzotov et al. I l|2004) . and was based on the Early Data 



Release. It studied the abundance patterns of metal-poor 
emission-line galaxies, with emphasis on the N/0 ratio. 
The second one is by iKniazev et al. I l(2fl04l) . who pro- 
duced a catalogue of emission-line galaxies with low oxy- 
gen abundances from the Data Release 1. In both cases, 
th e abundances were c omputed using the formulae given 
bv llzotov etlD 1)1994 . 

Here, we use the Data Release 3 ( DR3). We adopted se- 
lect ion criteria diffe r ent fr om those bv llzotov et al. I ( 2004}> 
and lKniazev et aL so that our sample is in general 

composed of objects with higher quality spectra and there- 
fore more reliable abundances. In the present paper we 
compute the abundances with revised expressions, based 
on the most recent atomic coefficients for the line emissiv- 
ities and on appropriate series of photoionization models 
for the ionization correction factors. All the results are 
available in electronic tables. We believe that they should 
provide a useful basis for many future works dealing with 
such objects. 

The organization of the paper is as follows. The selec- 
tion criteria used to extract the objects and the resulting 
sample are described in section |21 The derivation of the 
elemental abundances is presented in sectional The ob- 
tained abundance patterns are discussed in section^ The 
main conclusions of this study are presented in section [S| 



2. The samples 

2.1. The SDSS DR3 data 

The SDSS DR3 l)Abazaiian et al. 1 12005|) provides spec- 
tra in the wavelength range - 3800 - 9300A of 530 000 
galaxies, quasars and stars selected over 4188 square de- 
grees, and tables of measured parameters from these data. 
We have extracted the flux-calibrated spectra of the SDSS 
DR3 galaxies from the SDSS web page with the address 
^http: / /www . sdss.org/dr3 

From a visual examination of all the spectra and subse- 
quent measurements of emission line fluxes, we extracted 
~ 2700 spectra with an [O ill] A4363 emission line detected 
at a level better than la. High-excitation Sy2 galaxies 
with strong [O ill] A4363 emission line can be recognized 
because of their strong [O ii]A3727, He 11 A4686, [O ill] 
A5007 and [N 11] A6583 emission lines, and have been re- 
moved from the sample. The reduction of the spectra and 
the measurement of the line intensities is done following 
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Fig. 1. Diagnostic diagram [O ill] A5007/H/3 vs [N 11] 
A6583/Ha for our SDSS sample (See Sect 2.2). The dashed 
line separates galaxies with H 11 regions-like spectra from 
galaxies containing an a ctive nucleus (labeled "AGN") 
(iKauffmann et al. L l2003f) . As in all the figures with ob- 
servational data, the number of points is indicated in the 
left lower corner. 
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Fig. 2. Difference between ionic abundances calculated 
with the new and old emissivities in the case of the low- 
density limit as a function of the electron temperature. 
Curves are labeled by the emission lines used for the abun- 
dance determination. 



the prescriptions of llzotov et al. 1 1)2004|) . In particular, us- 
ing the IRAF/SPLOT routine, the line intensities are fit- 
ted with a gaussian profile, and the continuum is visually 
defined on both sides of each line. 

In the case of hydrogen lines, where the nebular emis- 
sion line sits on top of a much broader stellar absorption 
line, the level of the continuum is set at the bottom of 
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Fig. 3. Relations between the electron temperatures (a) tc(0 ii) and tc{0 ill) (tp = 10 ^ Tp) and (b) tpfS ill) and to(0 



III) obtained from the se quences of ph otoionization models for H ii galaxies ( Stasinska fc Izotovl 1200,'^ using stellar 
atmosphere models from lSmith et alH (|20021 . Filled circles are for the models with Z = 0.02 Zq, open circles are for 
the models with Z = 0.05 Zp, an d stars are for the models with Z ~ 0.2 Zq. Th e dotted l i nes ar e the approximation 
of the models of IStasihska I l)l990f) in (a) and the approximation of the models of iGarnettl l)l992j) in (b). 
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Fig. 4. (a) Relation between the electron temperatures tc{0 ll) and tc{0 ill) for the SDSS sample (small circles) and 
HeBCD galaxies (large circles), (b) Relation between the electron temperatures tc{S ill) and tc{0 ill) for the SDSS 
sample. The solid and dashed lines in (a) and (b) are respectively the predictions for H ii region models fEas. HM- 
I15|) with low- and intermediate metallicities. 



the absorption feature. Given the spectral resolution of 
the SDSS, this absorption feature is clearly seen whenever 
it is important. We estimate that, for the objects studied 
in this paper, which have large equivalent widths of the 
hydrogen emission lines, the presence of this underlying 
absorption introduces only a small uncertainty (< 10%) 
on the fluxes of the strong hydrogen emission lines. 

The 3800 - 9300A wavelength coverage of the SDSS 
spectra introduces two sources of uncertainties in the emis- 
sion line fluxes. First, the emission line [O ii] A3727 in the 
spectra of galaxies with z ^ 0.02 is not observed. The 
fraction of such galaxies is ^ 30%. In those cases the 
flux of the [O ii] A3727 emission line is estimated from 
the fluxes of the [O ii] A7320, 7330 emission lines tak- 



ing into account the empirical relation between the [O 
il] and [O III] temperatures (see Eq. I14II . iKniazev et al. I 
('2004') have shown that the use of the [O ii] A7320, 7330 
emission lines instead of the [O ii] A3727 emission does 
not introduce biases in the derived oxygen abundances. 
Therefore we use the former lines to derive the O"*" abun- 
dance in the H ii regions with the non-detected [O ii] 
A3727 emission line. Second, differential refraction effects 
may be important. The majority of spectra in our sample 
have been obtained at airmasses ^ 1.3. For an airmass of 
1.3, the offset of the H ii region (or of a standard star) 
at wavelengths 4000A and 7500A relative to its nominal 
location in side the slit at 50 00A are -|-0'.'5 and — 0'.'5, re- 
spectively llFilliDenkoLll982l) . For comparison, the radius 
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of the round slit used in the spectroscopic observations is 
1'.'5. Therefore, the observed fluxes of the blue lines [O ii] 
A3727 and [Ne ill] A3868 and the red Hues [Ar ill] A7135 
and [O ii] A7320, 7330 used for abundance determination 
could be affected by differential refraction. From existing 
data we cannot estimate the uncertainties introduced by 
this effect. Fortunately, strong hydrogen emission lines in 
the wavelength range ^ 3700 - 66OOA are present in the 
spectra of the majority of the selected galaxies. Correction 
for extinction using the observed fluxes of these lines also 
corrects at the same time for the effects of differential re- 
fraction. 

The precision of abundance determinations is critically 
dependent on the uncertainties in the measurements of 
the [O III] A4363 emission line fluxes. In the case of galax- 
ies with weak [O iii] A4363 e mission, observational uncer- 
tainties ijlzotov et al. may create spurious trends 
in the abundance patterns. Therefore, we restricted our 
initial sample of 2700 objects to objects with the best 
derived abundances. Empirically, we find that we can sig- 
nificantly reduce unphysical trends when restricting the 
sample to objects with an observed fiux in the H/3 emis- 
sion line larger than 10""'^'* erg s~'^cm~^. Additionally, we 
have excluded all galaxies with both [O iii]A4959/H/3 < 
0.7 and [O ii]A3727/H/3 > 1.0. A weak [O iii] A4363 line in 
these galaxies may be subject to significant enhancement 
due to a mechanism different from phot oionization heat- 
ing in H II regions ijlzotov et al. L l2004l) ^ . Applying all 
these selection criteria gives us a sample of ^ 310 SDSS 
objects, consisting mainly of high-excitation H 11 regions. 
The [O III] A4363 emission line intensities in these galaxies 
is detected at the 2a level for all of them and at the ^ 
5a for half of them. 

The observed line fiuxes F{X) are measured using the 
IRAF/SPLOT routine. The errors in the line fiuxes are 
calculated from the photon statistics. For this we use files 
attached to each spectrum which are generated by the 
pipeline data reduction and which contain the error in 
each pixel. Additionally, we adopt an error of 2% in the 
emission line fiuxes to account for the uncertainties of 
standard st ar absolute fiuxes used for fiux calibration (e.g., 
IOkeLll990t) . These errors will be later propagated into the 
calculation of abundance errors. The line fluxes were cor- 
rected for both reddening ( Whitford , 1958) and underly- 
ing hydrogen stellar absorption deriv ed simultaneously b y 
an iterative procedure as described in llzotov et al. Iljl994l) . 
This procedure automatically corrects emission line fluxes 
for differential refraction effects. The spectrum plate num- 
ber in the SDSS database, the extinction-corrected emis- 
sion line fiuxes /(A) normalized to /(H/3), together with 
their uncertainties, the equivalent width EW, the observed 



^ The second condition [O ll]A3727/H/3 > 1.0 ensures that, 
by removing objects with weak [O lll]A4959 we do not remove 
at the same time the most metal-deficient high-excitation H 11 
regions. For example, both the SE and NW components of I Zw 
18 have [O ni]A4959/H/3 < 0.7, but also weak [O n]A3727/H/3. 
Therefore, they are not removed from the sample. 
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Fig. 5. Ionization correction factor /CF(S+-1-S^+) vs the 
fraction 0+/{0++0'^+) for the "low models with dif- 
ferent rates of dielectronic recombination for sulfur ions. 
Nominal dielectronic recombination for sulfur ions is the 
same as that for oxygen ions. 



fiux F of the H/3 emission line and the extinction coeffi- 
cient C(H/3) are shown in Tabled 

Figurenshows the location of the galaxies of the SDSS 
sample in the [O ill] A5007/H/3 vs [N 11] A6583/Ha dia- 
gram which is generally used to distinguish H 11 galax- 
ies ionized only by massive stars from galaxies containing 
an active nucleus (AGN ). The dashed line separ ates the 
two classes of galaxies ( Kauffmann et ahl l2003l) . As in 
lizotov et al. (200^, all but one of our selected objects 
fall in the "H II" region of the diagram. 



2.2. The HeBCD sample 

In addition to the SDSS sample, we also use a sample of 
galaxies which has been collected primarily to study the 
helium abundances in low-metallicity blue compact dwarf 
g alaxies (the HeBCD sample). This sample is the same as 
in Izotov et al. '^('2004'), with the addition of galaxies from 
■Izotov & Thua,n . (..2004a.). 



3. Element abundance determination 

The physical conditions and element abundances are de- 
rived from the reddening-corrected emission line fiuxes 
using the electron temperature derived from the [O ill] 
A4363/(A4959-|-A5007) ratio. Temperature corrections to 
derive ionic abundances and ionization corrections to de- 
rive elemental abundances are based on sequences of pho- 
toionization models that best rep roduce the observed 
emission line trends in H 11 galaxies ijStasiiiska fc Izotov L 
i2003i) . 
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Table 1. Emission line fluxes in the SDSS galaxies^ 



Spectrum 


[Neiii] 


[Oiiil 


[Feiiil 

L J 


[Oiiil 


He I 


[0 il 


fSiiil 

L J 


[Niil 

L J 


[Sill 

L J 


[Arm] 


EW 




C(H/3) 


number 


3869 


4363 


4658 


4959 


5876 


6300 


6312 


6584 


6725 


7135 


A 






0266-51630-100^ 


0.370 


0.116 


0.005 


1.911 


0.117 


0.024 


0.017 


0.036 


0.142 


0.055 


262 


502 


0.00 




0.012 


0.005 


0.002 


0.057 


0.004 


0.002 


0.001 


0.002 


0.004 


0.003 








0267-51608-384* 


0.444 


0.042 


0.010 


1.263 


0.088 


0.020 


0.010 


0.096 


0.219 


0.053 


57 


314 


0.19 




0.018 


0.006 


0.004 


0.042 


0.005 


0.003 


0.003 


0.005 


0.008 


0.003 








0270-51909-306* 


0.175 


0.015 




0.967 


0.130 


0.026 


0.011 


0.294 


0.518 




55 


104 


0.10 




0.014 


0.009 




0.037 


0.010 


0.006 


0.006 


0.015 


0.019 










0272-51941-365* 


0.297 


0.023 


0.006 


1.131 


0.112 


0.058 


0.023 


0.224 


0.502 


0.076 


31 


110 


0.13 




0.019 


0.012 


0.011 


0.042 


0.010 


0.009 


0.008 


0.013 


0.020 


0.008 








0275-51910-429* 


0.430 


0.094 




1.587 


0.112 


0.029 


0.012 


0.059 


0.237 


0.056 


107 


197 


0.03 




0.017 


0.007 




0.051 


0.006 


0.004 


0.003 


0.004 


0.009 


0.004 









^Only the first five galaxies are shown here. The total list of galaxies 
He II A4686, [Ar iv] A4740, [Fc ill] A4988, [CI iii] A5517,5531, [O ii] 
here for lack of space, but they are present in the electronic table. 

■^in units 10^^^' erg s^^ cm^"^. 

*The first row are values, the second row arc errors. 



is available in electronic form. The hydrogen, [O II] A3727, 
A7320,7330 and [S iii] A9069 emission lines are not shown 
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Fig. 6. left Sulfur- to- oxygen abundance ratio vs the fraction 0+/(0++0^+) in the H ii regions with 12 -j- log 0/H 
< 7.6 for different rates of dielectronic recombination for sulfur. The nominal dielectronic recombination for S ions is 
the same as that for O ions. Large filled circles are data for the HeBCD sample, dots are for the SDSS sample. Solid 
horizontal line shows the mean S/0 abundance ratio from the data with the nominal rate of dielectronic recombination 
for sulfur ions, middle The same as in the left panel but for H ii regions with 7.6 < 12 -|- log 0/H < 8.2. right The 
same as in the left panel but for H ii regions with 12 -|- log 0/H > 8.2. The total number of galaxies in each panel is 
swown in the lower left corner. 
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Fig. 7. The same as in Fig. |51 but for argon. Nominal dielectronic recombination for Ar ions is the same as that for 
Ne ions. 



3.1. Atomic data and ionic abundances 

In previous studies, we have used emissivities for emis- 
sion lines provided by lAller I l)l984|) . Since that time the 
collisional strengths for many transitions have been re- 
evaluated. He re, we use t he at omic data from the refer- 
ences listed in lStasinska 1 1)2005|) . 

The electron temperatures derived from the [O 
iii]A4363/A5007 line ratio using old and new emissivities 
differ by less than 1% in the temperature range 5000- 
20000K. Therefore, we derive the electron tempe rature 
from a n iterative procedure, using the equation from lAller I 



1.432 



log[(A4959 + A5007)/A4363] - log Ct ' 
where t = 10-''Te(O ill), 



(1) 



Ct = (8.44 - 1.09t + 0.5^^ - O.OSt^)!^^^, (2) 

and x=10-'^ Net-° ''. The density iVe(S ii) is derived by 
interpolation of the [S ii] A6717/A6731 hne ratio vs Ne{S 
ii) relation) with th e use of the collisional strengths from 
lR,amsbottom et aTl ((1996.^ . The derived is always 



smaller than 10'^ cm so that the term containing x is 
never important. 

In Fig. 121 we show the differences between the ionic 
abundances calculated with new and old atomic data as 
a function of the electron temperature. Each curve is la- 
beled with the emission line used for the ionic abundance 
determination. We use [N ii] A6548, 6584 for the determi- 
nation of the nitrogen abundance, [O ii] A3727 (or [O ll] 
A7320,7331 when [O ii] A3727 is not available) and [O ill] 
A4959, 5007 for the oxygen abundance, [Ne ill] A3868 for 
the neon abundance, [S ill] A6312 and [S ii] A6717, 6731 
for the sulfur abundance when detected, [CI ill] A5517, 
5531 for the chlorine abundance, [Ar ill] A7135 (and [Ar 
iv] A4740 if seen) for the argon abundance, and [Fe ill] 
A4658 for the iron abundance. In general, the new atomic 
data change the ionic abundances by not more than 0.06 
dex. The only exception is Ar'^+ where the new determi- 
nations a re more than 0.1 dex lower than the old ones. 
Following 'Stasinska T (2005) we do not attempt to ana- 
lyze the c ollisional strengths f or Fe^"*" and use the same 
ones as in lThuan et al. I l)l995j) . We obtain linear fits for 
the dependences shown in Fig. El and use these fits in the 
abundance determinations. The ionic abundances are then 
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derived from the equations: 

12 + log 0+/H+ = log + 5.961 + ^ - 0.40 logi 



12-KlogO+/H+ 



-0.034i + log(l + 1.35a;) 

log ^7320+A7330 ^ g gg^ 



-t- 



H/3 

-0.483 log i 



-O.OlSt -I- log(l - 3.48a;), 



log A4959+A5007 ^ q 200 + 

-0.551og<-0.014t, 



12-hlog02+/H^ 

12-Hl0gN+/H+ = log ^6548+A6584 ^ 234 

_I_0J50 _ 0.42 log i 
-0.027t-Mog(l + 0.116a;), 

12 + log Ne^+/H+ = log + 6.444 + 

-0.42 logi- 0.009f. 



12-MogS+/H^ 



log ^6717+A6731 5 



0.929 



0.28 log t 



-0.018i-Mog(l-M.39x), 



(3) 



12-f logS'+/H+ = 

12 + logCl^+/H+ : 

(4) 12-hlogAr2+/H+ 
(5) 



log 2^ + 6.690 -t- 
-0.471ogt-0.010t, 



A5517-I-A5537 



5.700 

0.4831ogt- 0.002i 



1.088 



log 



A7135 



6.174 



H/3 

+ _ 0.48 logt 
-0.013t-f-log(l-h 0.22a;), 



12 + log Ar3+/H+ = log 2^ + 6.308 



1.280 



(6) 
(7) 

(8) 



H/3 ' ' t 

0.48 logt-0.047i. 



12-|-logFe^+/H+ = log 



A4658 
H/3 



6.498 4 
-0.48 logt. 



1.298 



(9) 



(10) 



(11) 



(12) 



(13) 



Since different ions reside in different parts of the H 
II region, the electron temperatures in Eq. |21-^]are not 
the same for the various ions. Only To(0 ill) is derived di- 
rectly from observations for all galaxies. Other-wise, we 
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Fig. 9. Ionization correction factors for heavy elements. Filled circles are for the models with Z — 0.02 Zq^ open 
circles are for the models with Z = 0.05 Zc:^. an d stars are for the models with Z = 0.2 Zq. The dotted line is the 
approximation of the models of IStasihska I l|l990t) . 



use the relations between To(0 ill) and the tempera- 
tures characteristic of other ions in sequences of pho- 
toionization models that best fit the observed line emis- 
sion t rends in H ii galaxies. The se sequences are defined 
as in IStasinska fc Izotov I l(200,'^ . but have been recom- 
puted with the new atomic data, and wi th an input radi- 
ation field computed with Starburst 99 l)Leitherer et al. I 
19991) using t he st ellar model atmospheres described in 
Smith et al. We thus obtain the relations between 

the electron temperatures Te(0 ill) and Te(0 ii) which are 
shown by solid lines in Fig. Eb- A dopting the nomencla- 
ture of IStasihska &: Izotov I l|2003l) , "low Z" models are 
represented by filled circles, "intermediate Z" models by 
open circles and "high Z" models by stars. These Z cor- 
respond to oxygen abundances 12 + log 0/H = 7.2, 7.6 
and 8.2, respectively. For comparis on, the dotted hne is 
t he fit of llzotov et al. I l)l994ll997a|) to the earlier models 
of IStasinskaTijiggOl) . The latter approximation, used in all 
our previous studies, is not very different from the new 
relations shown by solid lines in Fig. |3 We approximate 



the relations from the new models by the expressions 

t{0 II) = -0.577 + t X (2.065 - 0.498t), low Z, 

= -0.744 + t X (2.338 - 0.610i), intermed. Z, 
= 2.967 + t X (-4.797 + 2.827t), high Z, (14) 

where t{0 ii) = 10-'^T(O ii) and t = 10-'^T(O ill). 

In the H ii regions where both [O ii] A3727 and [O ii] 
A7320, 7331 are detected, the electron temperature To(0 
ii) can be derived directly from the A3727/(A7320+A7331) 
ratio. In Fig.0^ we compare the relations between Tc(0 ii) 
and rc(0 III) obtained from the observations with those 
from model predictions. Large filled circles correspond to 
galaxies from the HeBCD sample, while dots correspond 
to SDSS galaxies. Despite a large scatter of the data points 
caused mainly by large fiux errors of the weak [O ii] A7320, 
7331 emission lines, the relation between To(0 ii) and 
Te(0 III) derived from the observations follows generally 
the one obtained from the models. Therefore, in the fol- 
lowing we adopt the model relation to derive Tc(0 ii). 
Note that the scatter in Fig. ^implies that, for the low 
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redshift objects of the SDSS sample for which [O ii] A3727 
is outside the observed wavelength range, the 0+ abun- 
dance is not very accurate since it is obtained from [O ii] 
A7320,7331. This concerns about 30% of the sample. 

Similar relations between Tc(S ill) and Te(0 III) can 
be obtained from the same photoionization models: 

t{S III) = -1.085 + t X (2.320 - 0.420t), low Z, 

= -1.276 + t X (2.645 - 0.5m), intermed. Z, 
= 1.653 + t X (-2.261 + 1.6G5i), high Z, (15) 

where t{S ill) = 10-'*re(S ill) and t = 10-'*re(O ill). They 
are shown in Fig. ^p. The large wavelength coverage of 
the SDSS spectra allows us to detect at the same time the 
auroral [S ill] A6312 and nebular [S ill] A9069 emission 
lines in low redshift galaxies and to derive the electron 
temperature rc(S ill). In Fig.0lD we show the dependence 
of To(S III) on To(0 III) for those galaxies. Again, the 
theoretical relation is in satisfactory agreement with the 
observations, if we take into account the large scatter in 
the data, mostly due to measurement errors of the weak 
[S III] A6312 line. 

We thus adopt Te{0 ii) from Eq.^^for the calculation 
of N+, 0+, S+ and Fe2+ abundances, re(S ill) from Ea.[T31 
for the calculation of S^"*", Cl^+ and Ar^+ abundances and 
Te(0 III) from Eq.Qfor the calculation of 0^+ and Ne'^+ 
abundances. We adopt the "low Z" fits for H ii regions 
with 12 + log 0/H < 7.2, the "high Z" fits for the H ii 
regions with 12 + log 0/H > 8.2 (Eqs. [H and [13), and 
linearly interpolate Te{0 ll) and re(S ill) between "low Z" 
and "intermediate Z" cases for H ii regions with 12-|-log 
0/H between 7.2 and 7.6, and between "intermediate Z" 
and "high Z" cases for H ii regions with 12-1- log 0/H 
between 7.6 and 8.2. Since the population of the Ar'^+ ion 
is generally smaller than that of the Ar^+ ion and the [Ar 
iv] lines weaker than the [Ar ill] lines, we decided to derive 
the Ar abundance only from the abundance of the Ar^+ 
ion. 



3.2. Ionization correction factors 

The total element abundance can be derived from the 
abundances of ions seen in the optical spectra using ion- 
ization correction factors (ICFs). For this, the ionization 
structure of the H ii region has to be derived from pho- 
toionization models. Fortunately, the most abundant ions 
of oxygen, 0+ and 0^+, are seen in the optical spectra 
of the H II regions. This fact allows one to immediately 
derive the oxygen abundance as 



O 
H 



0+ O^ 
H+ ^ H^' 



(16) 



The small fraction of the unseen 0'^+ ion in the high- 
excitation H II regions can be added to the oxygen abun- 
dance if the He ii A4686 emission line is detected in their 
spectra. Our sequences of photoionization models show 
that O^^/O is ^1% only in the highest-excitation H ii 



regions with 0+/(0+-|-0^+) ^ 0.1. In those H ii regions, 
0^+/H+ can be approximated by 

Hc^+ /0^ 
Hc+ + Hc2+ I H^ 



O 



H^ 



= 0.5 X 



02 



H^ 



(17) 



This results in a minimal correction to the oxygen abun- 
dance derived just from 0+ and 0^+. In the lower- 
excitation H II regions with O'^ /{O'^+O'^^) > 0.1, we set 
0^+/H+ to zero. 

The ionization structure of an H ii region is deter- 
mined by the balance between the processes of photoion- 
ization and recombination. The rates for some of these 
processes are poorly known. One of the major problems 
is the uncertain rate of the dielectronic recombination for 
sulfur, chlorine and argon ions. The comparison of the 
abundances of these elements calculated with different as- 
sumed dielectronic recombination rates allows us to put 
some constraints on these rates. Fig. [S] shows the depen- 
dence of /CF(S++S2+) with 0+/(0+-f02+) for the "low 
Z" models. The solid curve corresponds to what we call 
the "nominal" dielectronic recombination rates (equal to 
the dielectronic recombination rates for second-row ele- 
ments in the same isoelectronic sequence). We adopt as 
the nominal dielectronic recombination rates for S ions 
the same as for O ions. The dotted line corresponds to 
dielectronic recombination rates enhanced by a factor 3, 
and the dashed one is obtained with a total suppression of 
dielectronic recombination. It is seen from this figure that 
/CF(S+-|-S2+) decreases when the dielectronic recombi- 
nation rate increases and the effect is more important for 
high-excitation models. Figure El shows the values of S/0 
calculated with different dielectronic recombination rates 
as a function of the fraction 0+/(0+-f 0^+) for the galax- 
ies in the SDSS and HcBCD combined sample. Since the 
effect of dielectronic recombination may depend on the 
metallicity in the H ii region, via the electron tempera- 
ture, we split the sample into three metallicity bins: the 
hottest H II regions with 12 + log 0/H < 7.6 (left panel), 
the intermediate-temperature H ii regions with 7.6 < 12 -|- 
log 0/H < 8.2 (middle panel) and the coolest H ii regions 
with 12-1- log 0/H > 8.2 (right panel). It is expected that 
if the adopted rate of dielectronic recombination for sul- 
fur ions is correct then there should be no dependence of 
the derived S/0 with 0+ /{0++0'^+). Inspection of Fig.lHl 
shows that trends in the S/0 abundance ratio arc minimal 
if a zero rate of dielectronic recombination for sulfur ions 
is adopted for the low- and high-metallicity bins and the 
nominal rate is adopted for the intermediate-metallicity 
bin. Thus, we adopt the ICFs{S+ +S'^+) calculated with 
a zero rate for dielectronic recombination for the low- 
and high-metallicity bins and with a nominal rate for the 
intermediate-metallicity bin. 

In a similar way, we consider the effect of dielec- 
tronic recombination on the ICFs for argon. We adopt 
as the nominal dielectronic recombination rate the same 
values as for ions of neon. Figure |7| shows the depen- 
dence of the Ar/0 abundance ratio with 0^/(0^+0'^^). 
Similarly to S/0, the value of Ar/0 should not depend 
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Table 2. Element abundances in the SDSS galaxies"'^ 



Spectrum number 


te{0 III) 


12+log O/H 


logN/O 


logNe/O 


logS/0 


logAr/O 


logFe/O 


0266-51630-100 


1.5533±0.0327 


7.84±0.02 


-1.62±0.03 


-0.85±0.04 


-1.68±0.04 


-2.40±0.04 


-1.93±0.14 


0267-51608-384 


1.2055±0.0600 


7.99±0.05 


-1.30±0.07 


-0.55±0.10 


-1.71±0.11 


-2.43±0.07 


-1.69±0.21 


0270-51909-306 


0.9444±0.1699 


8.33±0.20 


-1.09±0.26 


-0.77±0.41 


-1.67±0.33 






0272-51941-365 


1.0133±0.1661 


8.38±0.18 


-1.40±0.23 


-0.65±0.36 


-1.60±0.29 


-2.50±0.22 


-2.32±0.80 


0275-51910-429 


1.5117±0.0513 


7.82±0.03 


-1.45±0.05 


-0.72±0.06 


-1.73±0.10 


-2.37±0.05 





^Only the first five galaxies are shown. The total list of galaxies is available in electronic form. The Cl/O abundance ratios are 
not shown here for lack of space, but they are present in the electronic table. 



on 0+/(0+-|-0^+). A minimal trend in the Ar/0 abun- 
dance ratio is seen for the H ii regions in the whole range 
of metallicities if the nominal rate is adopted. A simi- 
lar consideration for chlorine (Fig. IS)) shows that minimal 
trends are present for the H ii regions with 12-1- log O/H 
> 7.6 if a zero rate is adopted. As for lower-metallicity 
H II regions, the statistics are too poor to make definite 
conclusions. We adopt above as the nominal dielectronic 
recombination rates for CI ions the same as for O ions. In 
the following, we adopt the ICFs derived from the mod- 
els with a nominal rate for dielectronic recombination for 
argon ions and a zero rate for chlorine ions. 

The adopted ICFs to obtain the elemental abundances 
from the ionic ones are shown in Fig. [SJ as a function of 
0+/(0++02+) or02+/(0++02+). AsinFig.El filled cir- 
cles, open circles and stars correspond to the H ii region 
models with Z = 0.02 Zq, 0.05 Zq and 0.2 Z©, respec- 
tively. Dot t ed lin es are approximations of ICFs from the 
IStasinska 1 l)l990l) models. 

It is seen from Fig. 1^1 that there is very little 
dependence of ICF(N+), ICF{Cl^+), ICF{kr^+) and 
/CF(Fe^"'") on metallicity over the whole range of oxy- 
gen abundance considered here. On the other hand, 
/CF(Ne2+) and ICF{S++S^+) are different for the Z = 
0.02 and 0.05 Zq model sequences as compared to the Z 
= 0.2 Zq sequence. This is due to the fact that the Z = 
0.02 and 0.05 Zq model sequences include an X-ray ioniz- 
ing source in addition to the ionizing star population. This 
was required in order to account for the large proportion 
of galaxies with nebular He ii A4686. If an X-ray source 
is present, it creates in the H ii region an extended warm 
zone where hydrogen is partly recombined, and where O^^ 
ions are decoupled from Ne^+ ions due to the efficiency of 
the 0^+ -I- H° charge transfer. At present, we have no real 
clue on the nature and strength of these X-ray sources, so 
the ICF for neon in low excitation H II regions with 12 + 
log O/H < 8.2 is uncertain. Fortunately, such cases do not 
represent a large fraction of our sample, as can be seen in 

Figs. El- El 

The new ICFs for nitrogen and iron are undisti n guish- 
able from the ones obtained from the IStasihskal l)l990l) 
models, while for other elements there are significant dif- 
ferences. We fit the new ICFs by the expressions 

/CF(N+) = -0.825U + 0.718 + 0.853/W, low Z, 

= -0.809u + 0.712 + 0.852/'i;, intermed. Z, 
= -lA76v + 1.752 + 0.688/i;, high Z, (18) 



ICF{Ne^+) = -0.385W + 1.365 + 0.022/w, low Z, 

= -0.405W + 1.382 + 0.021/u;, intermed. Z, 
= -0.591W + 0.927 + 0M6/W, high Z, (19) 

ICFiS+ + S2+) = 0.121U + 0.511 -t- 0.161/-!;, low Z, 
= 0.155W + 0.849 + 0.062/?;, inter. Z, 
= 0.178W -t- 0.610 -I- 0.153/w, high Z, 

(20) 

= 0.756W -t- 0.648 + 0.128/?;, low Z, 
= 0.814f; + 0.620 + 0.131/?;, intermed. Z, 

1.186?; + 0.357 + 0.121/?;, high Z, (21) 



/Ci^(Cl 



2+\ 



ICF{Ar'^+) = 0.278?; + 0.836 + 0.051/?;, low Z, 

^ 0.285?; + 0.833 + 0.051/?;, intermed. Z, 
= 0.517?; + 0.763 + 0.042/?;, high Z, (22) 



ICF{At^+ + Ar' 



0.158?; + 0.958 - 
= 0.104?; + 0.980 - 
0.238?; + 0.931 - 



0.004/?;, low Z, 
-0.001/?;, int.Z, 
0.004/?;, high Z, 
(23) 



/Ci^(Fe^+) = 0.036?; - 0.146 + 1.386/?;, low Z, 

= 0.301?; - 0.259 + 1.367/?;, intermed. Z, 
= -1.377?; + 1.606 + 1.045/?;, high Z, (24) 

where ?;=0+/(0++02+) and w=0^+ /{0++0^+). The 
fits given by Eqs. ^1 E| and |31 are applicable over the 
whole range of v,w — - 1, while the remaining fits are 
applicable only for v ^ 0.8. We note however that only 
a few H II regions in our sample have v > 0.8. We adopt 
the "low Z" ICFs for H ii regions with 12+log O/H < 
7.2, "high Z" ICFs for H ii regions with 12+log O/H > 
8.2 and linearly interpolated ICFs between "low Z" and 
"intermediate Z" cases for H ii regions with 12+log O/H 
between 7.2 and 7.6, and between "intermediate Z" and 
"high Z" cases for H ii regions with 12+log O/H between 
7.6 and 8.2. Although ICFs for some elements (e.g., for 
nitrogen) appear very similar in all metallicity bins (Fig. 
lUJ, we show in Eqs. fits obtained separately for 

each bin. 

In such a way, we can derive all elemental abun- 
dances using only analytical formulae. This is useful for 
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the derivation of the formal uncertainties in the derived 
abundances, since they can be obtained by expanding the 
formulae in Taylor series, and do not require the use of 
Monte-Carlo methods. 

4. Discussion 

4.1. SDSS DR3 galaxies versus galaxies from the 
HeBCD sample 

Electron temperatures te{0 ill), element abundances and 
abundance ratios in the H ii regions from the SDSS DR3 
sample calculated as explained above are shown in Table 
12 Some objects from the SDSS sample are also present 
in the HeBCD sample. For those objects we compare in 
Table 13 oxygen abundances derived from different obser- 
vations. The SDSS spectra of these galaxies are shown in 
Fig. ^1 Since the redshifts of these galaxies are small, 
the [O ii] A3727 is not in the observed wavelength range 
of the SDSS spectra and therefore the 0+ abundance 
was obtained from the [O ii] A7320, 7331 emission line 
fluxes. Exceptions are HS 0837-f 4717, HS 0924-f 3821, SBS 
0948-^532 and SBS 1358-f 576 with higher redshifts where 
the 0+ abundance was obtained from the [O ii] A3727 
emission line fluxes. The galaxies HS 0837-f4717, Mrk 
1236, SBS 0917+527, UM 461 and UM 462 were observed 
twice in the course of SDSS and we show both oxygen 
abundances for each of those galaxies. It is seen from Table 
Othat 12-|-log 0/H derived from different observations are 
generally in good agreement. One exception is the galaxy 
UM 461. In one of two SDSS observations of this galaxy, 
the [O III] A4959, 5007 emission lines are clipped resulting 
in a low oxygen abundance. The oxygen abundances for 
some galaxies from Table El are also in agreement with in- 
dependent abundance dete r minations fr om the same SDSS 
spectra bv lKniazev et al. I l)2003ll2004j) . 

In Fig. nil we show the computed N/0, Ne/0, S/0, 
Cl/0, Ar/0, Fe/0 abundance ratios vs oxygen abundance 
12 + log 0/H for the SDSS and HeBCD merged sample 
using old (left panels) and new (right panels) expressions 
to derive the abundances. In all panels, the HeBCD galax- 
ies are represented by large filled circles, while the SDSS 
galaxies are shown by dots. In each panel, t he solar abun- 
dance ratio (as compiled by lLoddersll2003^ is represented 
by a large open circle. 

It is seen from Fig. nil that SDSS galaxies are m gen- 
eral more metal-rich than the galaxies from the HeBCD 
sample. Only a few SDSS H ii regions have an oxygen 
abundance 12-|-log 0/H smaller than 7.6. Two of them are 
actually the northwest and the southeast components of 
the well-known metal-deficient star-forming dwarf galaxy 
I Zw 18, which are also present in the HeBCD sample. 
Thus, star- forming galaxies with extremely low metallic- 
ity are very rare in the spectroscopic data base of the 
SDSS. 

The distribution of the SDSS DR3 galaxies in the dif- 
ferent panels of Fig. | ll | con firms previous findings concern- 
ing HeBCD l|Thuan et al. . .1995: Izotov & Thuan ...1999,) 



Table 3. Comparison of element abundances in the galax- 
ies present in both HeBCD and SDSS samples 



Name (12+IogO/H)sDss (12-|-logQ/H)HcBCD 





7.84 


-|- 


0.02 


7.78 


-|- 


0.01 


J- J- kJ VJ O .i; 'J ^ 


7.42 


-|- 


0.03 


7.45 


-|- 


0.01 




7 (SS 


-j- 


n ns 


7 fil 


-|- 


m 




7.60 


-|- 


0.03 










8.15 


-|- 


0.04 


8.09 


-|- 


0.02 


T7w1SNW 

J.ZJVV-LOi'l VV 


7.14 


-|- 


0.05 


7.16 


-|- 


0.02 


T7wl8SF, 


7.25 


± 


0.04 


7.19 


-|- 


0.03 


IVli iVtJiJ 


8.29 




0.12 


8.37 


-|- 


0.01 


Mrki 7S 


7.90 


-|- 


0.11 


7.82 


-|- 


0.02 


Mrkl2Sfi 

i.V±l. XV X J-j O w 


7.99 


-|- 


0.05 


8.09 


-|- 


0.03 




8.12 


-|- 


0.02 








SRSnqi 7-1-^97 


7.81 


-|- 


0.04 


7.87 


-|- 


0.02 




7.91 


-|- 


0.02 








SRSnQ9fi4-finf! 


7.99 


-|- 


0.02 


7.92 


-|- 


0.02 




8.04 


-|- 


0.02 


8.01 


-|- 


0.01 


SRSnQ4S-l-5S2 

kJ -LJ kJ C t: LI "J 'J ^ 
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^ [O III] A4959, 5007 emission lines are clipped. 



and SDSS Early Data Release galaxies l|lzotov et al. L 
\2()04\ . The use of the new expressions to compute elemen- 
tal abundances (FigE] middle and right panels) result in 
very slight changes. Ne/0 abundance ratios are typically 
decreased by 0.1 dex. The largest effect is for the S/0 and 
Ar/0 abundance ratios, which decrease typically by 0.15 
dex. 

The mean Ne/0 and Cl/0 abundance ratios for our 
galaxies ar e in good a greein ent with the solar abundance 
ratios from iLodders I (ji^Q^) (open circles in Fig. [TT|l . On 
the other hand, the S/0 and Ar/0 abundance ratios in the 
H II regions are lower than the solar abundance ratios (by 
0.15 dex and 0.25 dex, respectively). We cannot exclude 
the possibility that these differences are caused by uncer- 
tainties in the determination of S and Ar abundances in 
the H II regions. Alternatively, they may indicate that the 
solar S and Ar abundances are not as accurately known 
as thought. 

We show by solid lines in Fig. ^2 (middle and right 
panels) the regression lines for log Ne/0, log S/0, log 
Cl/0, log Ar/0 and log Fe/0 vs 12-hlog 0/H obtained 
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Fig. 10. SDSS spectra of the H ii regions listed in TableO The names of the spectra in the SDSS database are shown 
in parentheses. 



from the HeBCD sample only. The equations of these lines -^og — — —(0 026 =b 027)X — (1 514 zb 215) (26) 
are the following: ^ 



Ne 

log — = (0.088 =b 0.018)X - (1.450 =b 0.144), 



O 

(25) log — - (0.120 =b 0.106)X - (4.365 =b 0.852), 



(27) 
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Fig. 11. log N/0 (a, g), log Ne/0 (b, h), log S/0 (c, i), log Cl/0 (d, j), log Ar/0 (e, k) and log Fe/0 (f, 1) vs oxygen 
abundance 12 + log 0/H for the emission-line galaxies. Large filled circles are galaxies from the HeBCD sample, 
dots are galaxies from the SPS S sample. Abundanc es in the left panel are calculated with the expressions used by 
llzotov et al. I () 19941 Il997al) and iThuan et al. I lll99fil) and th ose in the right panel are obtained with expressions from 
this paper. The solar values as compiled bv lLodders 1 1 200,'^ are indicated by the large open circles and the associated 
error bars are shown. Straight lines are linear regressions obtained for the HeBCD sample. 



log — = (0.039 ± 0.031)X - (2.714 ± 0.244), (28) 
Fe 

log — = -(0.606 ± 0.081)X -I- (2.994 ± 0.639), (29) 

where X = 12 + log 0/H. 

For the merged HeBCD and SDSS sample we find no 
significant trends with the oxygen abundance for the S/0, 
Cl/0 and Ar/0 abundance ratios in the whole range of 
oxygen abundance (Figs. Illh - lllk). The best derived 
abundances among the a-elements are for oxygen and 
neon because the atomic data for these elements are better 
known, the abundances are derived from bright emission 
lines and only two ionization stages are present in signifi- 
cant amount in H ii regions. 

4.2. Elemental depletion onto dust grains 

We note that there is a slight increase in Ne/0 with in- 
creasing 0/H in panel h of Fig. (the correlation is sig- 
nificant at the 99.98% level according to the Spearman 
statistical test). If the new expressions for deriving abun- 
dances are correct then this slight increase (by ^0.1 dex 
over the entire metallicity range considered) can be in- 
terpreted by depletion of oxygen onto dust grains. This 



corresponds to ^ 20% of oxygen locked in the dust grains 
in the highest-metallicity H ii regions of our sample, while 
no significant depletion would be present in the H ii re- 
gions with lower metallicity (assuming that the adopted 
solar Ne/0 ratio corresponds to that of the interstellar 
medium). Since Ar is a noble gas like Ne, one would ex- 
pect the same trend in Ar/0 as in Ne/0. The regression 
lines for the merged HeBCD and SDSS sample shown in 
Fig. Illh and lllk are 

Ne 

log — = (0.097 ± 0.015)X - (1.531 ± 0.120), (30) 
Ar 

log — = -(0.084 ± 0.021)X - (1.729 ± 0.166), (31) 

It is seen that the slopes of the Ne/0 vs 0/H and Ar/0 
vs 0/H regressions are statistically different. This could 
indicate that the production of Ar relative to O by the 
stellar populations decreases with increasing metallicity. 

As seen in Fig. 1111 iron has a very different behavior 
as compared with other elements. The Fe/0 ratio is well 
below solar for most of the galaxies. In addition, the trend 
with 0/ H is opposite to that in metal-poor stars in our 
Galaxy I' Carretta etini200dlTakedaLl20o3lBai et all 
i2004) . This argues in favor of depletion of iron onto dust 
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Fig. 12. (a)— (c) Nitrogen-to-oxygen abundance ratio, log N/0, vs equivalent width of H/3 emission line, EW(H/3), 
for the galaxies in three metallicity bins, (d)— (f) The same as in (a)— (c), but for neon-to-oxygen abundance ratio, 
(g)— (i) The same as in (a)— (c), but for sulfur-to-oxygen abundance ratio, (j)— (1) The same as in (a)— (c), but for 
argon-to-oxygen abundance ratio, (m)— (o) The same as in (a)— (c), but for iron-to-oxygen abundance ratio. Straight 
lines are linear regressions obtained for intermediate Z and high Z samples. 



grains. The decrease of the Fe/0 abundance ratio with 
increasing 0/H (by ~ 1 dex when the oxygen abundance 
increases by ~ 1.5 dex as seen in Fig.lllL^ implies that de- 
pletion of Fe increases with increasing metall i city, s imilar 
to what is suggested bv iRodriguez fc R.ubin I ll2005l) from 
the consideration of a small sample of Galactic and ex- 
tragalactic H ii regions. If this explanation is correct then 
only ^ 10% of iron in the high-metallicity H ii regions of 
our sample (log 0/H +12 — 8.2 - 8.5) is in gas phase and 
the remaining ^ 90% of iron is locked in dust grains, while 
in the H li regions with the lowest metallicity the depletion 
is much lower. This conclusion is supported by the small 
trend seen in the Ne/0 abundance ratio. Since Ne is not 
locked in dust grains, the increase of Ne/0 with metallic- 
ity can be explained by a small depletion of oxygen. Since 
oxygen is an element ^ 15 times more abundant than 
iron, a depletion of oxygen by ^ 20% would correspond to 
an important iron depletion. In the panels on the right of 
Fig. 1121 we show the variation of Fe/0 as a function of the 
equivalent width of H/3, EW(H/3), for the galaxies in the 
SDSS and HeBCD merged sample. Fig. I12h i corresponds 
to the lowest metallicity bin, defined by 12 -I- log 0/H < 
7.6, Fig. I12h corresponds to the intermediate metallicity 
bin (7.6 < 12 + log 0/H < 8.2) and Fig. IHb corresponds 
to the highest metallicity bin in our sample (12-1- log 0/H 
> 8.2). If we take E W(H/3) as a measure of the age of the 
ionizing stars (but see lStasinska fc Izotov llioosll . Fig.lT^ 
strongly suggests that the grains are gradually destroyed 
on a time scale of several Myr. The Spearman rank corre- 



lation coefficient is Rs— -0.30 with an associated signifi- 
cance level of 99.9%, indicating that the observed correla- 
tion is not due to chance. There is no similar clear trend 
in Fig. I12h {Rs= -0.006 with an associated significance 
level of 7%). However, since the initial depletion in the 
intermediate metallicity bin is lower than in the highest 
metallicity bin, detecting a temporal variation of the cor- 
relation is more difficult. Note that the scatter in Fe/H 
at any value of EW(H/3) as well as at any value of 0/H 
is far more important in Fe/0 than in Ne/0, S/0, Cl/0 
and Ar/0. We believe that this is mostly due to errors in 
the iron abundance determinations since the [Fe ill] lines 
are generally very weak but one cannot exclude that the 
degree of iron depletion might vary among H ii regions of 
similar characteristics. 

4.3. The nitrogen conundrum 

Much has already been written about the N/0 v ersus O /H 
diagram (upper line o f panels in Fig . HTll ( e.g.. iGarnett 1 . 
Il990t IPilvugin I Il993t iHenrv et all l2000|) . Our new re- 
sults basically confirm the findings in o ur previous papers , 
based on considerably sm aller samples ijizotov fc ThuanT 
ll999Hlzotov et al. LlioO^ : the N/0 ratio decreases as 0/H 
decreases but reaches a plateau at 12 -|- log 0/H < 7.6. 
Moreover, while at a given 0/H, there is a significant dis- 
persion in N/0 when 12 + log 0/H > 7.6, the dispersion is 
very small at lower values of 0/H. The lack of H ii regions 
with 12 + log 0/H < 7.6 and high N/0 is likely not due to 
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selection effects. Indeed higher N/0 abundance ratio im- 
phes higher flux of the [N ii] A6584 emission line which can 
be more easily detected in this case. On the other hand, 
it is unlikely that we missed low-metallicity galaxies with 
a very low N/0 abundance ratio, below the plateau value 
of log N/0 ~ -1.6, as the [N ii] A6584 emission line was 
detected in all most metal-deficient galaxies. 

The N/0 versus 0/H pattern for H ii regions differs 
signif icantlv from that for damped Lyg sy stems (DLA) 
(;e.e-.. ICenturi6n et al. ll20n.'4lMolaroll200,^ . in that some 
DLAs have N/0 ratios much lower than -1.6 dex, the 
value of the observed plateau for H ii regions. Recently, the 
N/0 - 0/ H diagram has also b een constructed for metal- 
poor stars (ISnite et al. L l200fil) . Stars affected by mixing 
have been excluded from their diagram. It is seen that 
some very low-metallicity stars also have values of N/0 
below the H ii region plateau, although they are larger 
than the lowest N/0 found for DLAs. To interpret these 
results in terms of models of chemical evolution of galax- 
ies, we should however bear in mind that abundances de- 
termined in DLAs fseelDessauges- Zavadsk v et al. LEo04I) 
as well as in stars ijSoite et al. 1 12005,,) are generally less 
accurate than the N/0 ratios in metal-poor H II regions 
where the main source of uncertainty is the measurement 
of the weak [N ii] 6584 hne. In the HeBCD sample, the 
average formal error in N/0 is ^ 0.1 dex at 12 + log 
0/H < 7.7, and it is < 0.2 dex for bright galaxies in the 
SDSS DR3 sample at the same metallicities. Also, DLAs, 
Galactic halo stars and metal-poor emission-line galaxies 
had likely very different star formation histories. The dif- 
ferences in the abundance patterns observed between these 
systems might also result from the different timescales and 
intensities of the star formation episodes which enriched 
their interstellar medium. 

The fact that no H ii region is found with an extremely 
low N/0 abundance ratio ( log N/0 < -1.6) suggests that 
the SDSS galaxies are not extremely young and have ages 
more than 100 - 300 Myr, which are required for the com- 
pletion of the evolution of the intermediate-mass stars. 

We also confirm the trend of increasi ng N/0 with 
decreasing EW(H/3) already mentioned by llzotov et al. I 
l|2004l) . This trend is seen in the 12 -f log 0/H > 8.2 
metallicity bin (at a 99.99% significance level according 
to the Spearman statistics) and also in the 7.6 < 12 -|- 
log 0/H < 8.2 metallicity bin (at a 99.8% significance 
level) (panels c and b of Fig. I12|l. Note that, apart from 
Fe/0, which we discussed before, none of the remaining 
abundance ratios shown in Fig. 1121 exhibits such a trend. 
Therefore, the trend for N/0 with EW(H/?) must be real 
and not due to selection effects or systematic errors in 
abundance determinations. 

Since EW(H/3) is a measure of the age of the ionizing- 
star population, it is interesting to check whether the 
larger N/0 abundance ratio in the more evolved starbursts 
is due to the nitrogen enrichment by the Wolf-Rayet stars. 

Let us consider the massive star rotating mod- 
els for a metallicity equal to Z=0.004, which cor- 
responds to 12-l-log 0/H=8.2, computed recently by 



iMevnet fc Maeder I l)2005|) . The 60Mq star model loses 
about half of its initial mass during the first 4.76 Myr. 
During this phase the amount of nitrogen in the wind 
ejected material is about 9 times the initial content (the 
initial content would be the mass of nitrogen ejected if 
no change of nitrogen abundance occurred) . On the other 
hand, oxygen appears to be depleted by a factor 2.5 with 
respect to its initial value. Thus, we conclude that if one 
were to observe only wind material of this kind then one 
would expect an increase of the N/0 ratio by about a fac- 
tor 20 in a few Myrs. A similar estimate made on the basis 
of a 120Mq stellar model leads to an N/0 increase by a 
factor 55, while a SOMq stellar model would predict an 
increase of only a factor 5. Of course, there is a dilution of 
the WR material with the local interstellar material of the 
H II region ^. The enhancement factor of the N/0 abun- 
dance ratio in the H ii region due to the WR wind can be 
estimated as 



N/0(Hn) _ p(N) X {M^/M,) X (M,/Af (HH)) + 1 
N/0(ini) ^ 39(0) X (M^/M,) x (M,/M(HH)) + r 



(32) 



where p(N) and p{0) are the production factors of nitro- 
gen and oxygen in the wind, is the total mass in the 
winds, Mt. is the integrated initial mass of the stars and 
M(H ii) is the mass of the H ii region. Let us take p(N)=9, 
p(O)=0.4, typical of the 6OM0 model. In the mass range 
between 30 and 120 M©, on average about half of the stel- 
lar mass is lost through stellar winds enriched in N and 
depleted in O, thus 



/3o'° 0.5AfM-2-35d7\f 
J^f MM-2-35dM 



0.03. 



(33) 



About 3 percent of the mass in stars is ejected under the 
form of N-rich stellar wind. For an H 11 region of constant 
number density n, the value of M*/M(H ll) for a clus- 
ter with a Salpeter initial mass function is roughly equal 
to n/360. If the number density of the H 11 region is 100 
cm^^, M»/M(H 11) is about 0.28. Thus the enhancement 
factor of the N/0 abundance ratio in the H 11 region is ^ 
1.07 (or ^ 0.03 dex), significantly lower than the disper- 
sion of the N/0 abundance ratio in Fig. 1111 and the value 
of the trend seen in Fig. ^] (left panel). 

It is likely, therefore, that the large dispersion of N/0 
abundance ratio in Fig. ^2 and the trend seen in Fig. E| 
are not due to a global increase of the N/0 ratio in the H 
II region, but it could still be a local effect. It was assumed 
in the above consideration that the H 11 region is uniform, 
i.e. its density is constant. However, it is very likely that 
the density of the N enriched ejecta is much higher than 
that in the ambient H 11 region. Since the luminosity of 
the forbidden lines in the low-density medium scales as 
the number density of the emitting ions multiplied by the 



^ Note also that these stars may also eject through their 
winds and at the time of the supernova explosion material with 
very low N/O ratios. This would also concur, together with the 
dilution effect, to decrease the N/O ratio in the H 11 region. 
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number density of electrons, the apparent enhancement of 
the N/0 abundance ratio, can be estimated from 



N/0(HII) 



N/0(ini) 



app 



p(N) X noj/n X M„/Af (HII) + 1 
p{0) X nej/n X M„/A//(HII) + 1 



,(34) 



where ricj is the number density of gas particles in the 



extremely low-metallicity galaxies as they were collected 
for the determination of the primordial He abundance. 
The small number of extremely low-metallicity galaxies in 
the SDSS may suggest that such galaxies are actually very 
rare. However, the selection effects in the SDSS might be 
important because most metal-poor galaxies are in general 
low-luminosity objects and the SDSS spectroscopic data 



ejecta. Adopting n = 100 cm and riej = 1000 cm we are onlv comp lete for galaxies brighter than ^ 17.7 mag 



obtain an apparent enhancement factor of ~ 1.7 from Eq, 
1341 comparable to the dispersion of the N/0 abundance 
ratio in Fig. II II and the value of the trend in Fig.^l Note 
that there is no direct diagnostic of the density of the 
emitting ejecta. The density derived from the [S ii] line 
ratio applies essentially to the H ii since the winds are not 
enhanced in sulfur. 

In summary, we find that the observed trend of N/0 
increasing as EW(H/3) decreases is naturally explained by 
the expected ejection from Wolf-Rayet stars. However, the 
true value of the N/0 enhancement factor is defined by 
Eq. 1221 and is small. The dispersion of the observed N/0 
abundance ratio is caused by the local enrichment by the 
WR stars and the value of the N/0 abundance ratio in the 
H II before enrichment by the present starburst is likely de- 
fined by the lower boundary in the N/0 vs 0/H diagram. 
If this is correct then the small dispersion of the N/0 
abundance ratio in the most-metal deficient H ii regions 
could perhaps be explained by the very low number of the 
WR stars and short WR stage during the low-metallicity 
burst of star formation. However, a larger sample of ex- 
tremely low-metallicity objects would be needed to con- 
firm this. 

5. Conclusions 

We have examined all the galaxies in the Sloan Digital 
Sky Survey (SDSS) Data Release 3 (DR3) to select out 
those with a detected [O ill] A4363 emission line, which 
allows, in principle, a direct element abundance determi- 
nation based on the electron temperature. To minimize 
errors in the abundance determinations, we have kept for 
further analysis only ^ 310 SDSS DR3 galaxies that are 
sufiiciently bright in H/3 (i^(H/3) > 10"" erg s"^ cm'^). 
The [O III] A4363 line in these objects is detected at the 
^ 2(7 level. This sample was merged with a sample of 109 
blue compact dwarf galaxies with high-quality observa- 
tions and large equivalent widths of emission lines, referred 
to as the HeBCD sample. In this way we obtained samples 
of emission-line galaxies, spanning a range in metallicities 
from 12 -I- log 0/H - 7.1 (Zq/30) to - 8.5 (0.7 Zq). Our 
main conclusions are as follows: 

1. Despite an examination of the entire SDSS DR3 
sample of ~ 530 000 galaxies, we found only 6 new galax- 
ies with extremely low metallicity (12 + log 0/H < 7.6, 
i.e. Z < Zf7>/12 ), in addition to the galaxies discovered by 
iKniazev et all (2003). For comparison, the HeBCD sam- 
ple which consists of only '--^100 galaxies contains about 15 
such galaxies. This difference is in part due to selection ef- 
fects. The selection criteria for the HeBCD sample favour 



are only complete to i 
l|Tzotov et a].Ll2004|>. 



2. The emission-line galaxies in the SDSS sam- 
ple show distributions of the N/0, Ne/0, S/0, Cl/0, 
Ar/0, Fe/0 abundance ratios vs oxygen abundance 
similar to those found earlier for the galaxies in the 
HeBCD sample llThuan et al. L Il995l llzotov fc ThuanL 
ll999HCxuseva et al. Ll2003[r 

3. The a element-to-oxygen abundance ratios Ne/0, 
Cl/0, S/0 and Ar/0 do not show large trends with oxy- 
gen abundance. The best determined ratio, Ne/0, in- 
creases slightly with increasing 0/H (by ~ 0.1 dex in the 
considered metallicity range). This can be explained by 
the depletion of oxygen onto dust grains with ~ 20% of O 
locked in dust. No significant depletion of oxygen is found 
in the lower-metallicity H ii regions. The mean Ne/0 and 
Cl/0 abundance ratios in the H ii regions from our sam- 
ple are consistent with the solar ratios. On the other hand, 
the mean S/0 and Ar/0 abundance ratios in the H ii re- 
gions are lower than the solar ratio. This may indicate the 
S and Ar abundances in the Sun are somewhat uncertain. 

4. The Fe/0 abundance ratio shows an underabun- 
dance of iron relative to oxygen as compared to solar, 
suggesting depletion of iron onto dust grains. In the high- 
metallicity galaxies more than 90% of iron is locked onto 
dust grains. On the other hand, no significant iron deple- 
tion is present in the lowest-metallicity galaxies. We also 
find evidence that dust destruction occurs on a time scale 
of several Myr. 

5. No galaxy with log N/0 ^ -1.6 was found, suggest- 
ing that the chemical evolution of local low-metallicity 
emission-line galaxies is different from that of some high- 
redshift DLAs which have considerably lower log N/0 (~ 
-2.3). These DLAs are considered to be truly young galax- 
ies, with nitrogen produced only by massive stars. If this 
interpretation is correct, then our sample of ~ 400 dwarf 
emission-line galaxies contains no extremely young galaxy 
but rather galaxies of ages ^ 100 - 300 Myr, required for 
the enrichment in nitrogen by intermediate-mass stars. 

6. Our data indicate an apparent increase of N/0 with 
decreasing EW(H/3), best seen among galaxies from our 
sample that have intermediate metallicities. We interpret 
this as evidence for gradual enrichment of the H ii region in 
nitrogen by massive stars from the most recent starburst. 
The magnitude of the observed effect is consistent with 
that expected from current theoretical models for massive 
stars, if account is taken for the density of the ejecta being 
larger than that of the H ii regions. 
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